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Abstract 

We discuss effects of the stepping kinetics of molecular motors on their traffic behavior on crowded 
filaments using a simple two-state chemomechanical cycle. While the general traffic behavior is 
quite robust with respect to the detailed kinetics of the step, a few observable parameters exhibit a 
strong dependence on these parameters. Most strikingly, the effective unbinding rate of the motors 
may both increase and decrease with increasing traffic density, depending on the details of the 
motor step. Likewise the run length either exhibits a strong decrease or almost no dependence on 
the traffic density. We discuss recent experiments in light of this analysis. 
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I. INTRODUCTION 



Molecular motors power the transport of various kinds of cargoes within cells by directed 
steppmg movements along filaments of the cytosUeleton (11]. The similarities (and differenees) 

compared to highway traffic [2| and the fact that cells are crowded environments l3|, 4 1 have 



stimu 



10 



ated extensive research of traffic phenomena in molecular motor systems [5|, l6 
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Most studies of molecular motor traffic have modeled the stochastic stepping of an indi- 
vidual motor by a single Poissonian step, which is modified in dense traffic by an exclusion 



rule similar to the well-studied asymmetric simple exclusion process (ASEP) |14j . Il5l . 
If the site to which a motor attempts to step is not accessible, because another motor (or any 
other kind of obstacle) is bound there, the step is rejected. Stepping of a molecular motor 
is however a complex process and consists of a series of transitions between different motor 
states and corresponding conformational changes |18|]. Several recent studies have therefore 
incorporated detailed kinetic models of the steps into models for molecular motor traffic 
201 1 ■ However, the main feature of molecular motor traffic visible in experiments 
22|, the emergence of traffic jam-like density profiles, with a region of high density 
of bound motors separated from a low-density region by a sharp interface 
as the basic features of the phase diagrams for transport in open systems [6|, 
robust with respect to such extensions of the simplest models. 

In this paper, we discuss observable effects of the stepping kinetics on the traffic behavior. 
To be explicit we consider the simplest possible stepping cycle, which consists of transitions 
between two motor conformations or internal states of the motor, with only one transition 
involving movement of the motor to the next binding site on the filament as shown in Fig. [H 
A few generalizations for models that incorporate reverse transitions and backward steps 
are briefly discussed in the appendix. 





23[ |, as well 

are very 



II. MODEL WITHOUT MOTOR BINDING/UNBINDING 

We flrst consider the case without exchange of motors with a reservoir, i.e. we neglect 
that fact that bound motors unbind from their fllamentous track and unbound motors bind 
to it. In the simplest model, the motor cycles between two conformational states, state 1 



2 



and 2, which may be considered as a weakly bound and a strongly bound state as in a kinetic 



scheme for kinesin 1 



24| . In that case, the model with two substeps is defined by two rates, 



Ki and K2. The first one describes the transitions from state 1 to state 2 with the motor 
position along the filament unchanged, while the second describes the actual step to the next 
site, which is accompanied by a transition from motor state 2 (the strongly bound state of 
kinesin) to state 1 (the weakly bound state) and has step size as depicted in Fig. [T](a). 
Binding and unbinding of motors as indicated in Fig. [^b) will be incorporated into the 
model in the next section. Throughout this paper we use periodic boundary conditions; we 
note however that more realistic open boundaries will lead to the same results if the system 
size is long compared to the run length of the motors, i.e. the distance a motor moves along 



the filament before unbinding 



2|]. Under the latter condition, which is typically fulfilled for 



cytoskeletal motors with run lengths of ~ 1/im and filaments with typical lengths of tens of 
microns the choice of the boundary conditions is only refiected in boundary layers of a 
size comparable to the run length Q|. 

Fig. El^a) shows the motor current as a function of the motor density p on the filament as 
obtained from simulations. Here and in the following plots, we normalize the motor current 
J and the motor velocity v by the single motor velocity Vq = k,iK,2^/{i^i + 1^2), which we 
consider as a known quantity given by experimental data. In contrast to the quantities 
J, V and Vq, the normalized current Ji/v^ and the normalized velocity v/vq, which will 
be discussed below, depend only on the ratio ki/k2 of the two transition rates and not 
on the absolute values of ki and k,2- Therefore JC/vq was plotted for different choices of 
K\j K2 in Fig. [2](a). In the limit k\ ^ /t2, the actual step is rate-limiting and the transition 
from state 1 to state 2, which does not comprise movement along the filament, can be 
neglected. In this limit, the two-step model becomes equivalent to the usual ASEP and 
the current is given by J = (fo/^)p(l — p), symmetric around p = 1/2. If however the 
conformational transition cannot be neglected, this symmetry is lost and the maximum of 
the current is shifted towards higher densities. At the same time, the value of the maximal 
current increases. The absence of particle-hole symmetry, i.e. symmetry upon interchanging 
p and (1 — p), can be understood by considering the velocity of rare individual particles or 
holes, i.e. the limits p ~ and p ~ 1, respectively. While a single motor particle moves with 
velocity fo, a single hole in a fully occupied lattice moves with the hole velocity vu = ^2^, 
as all particles it encounters had sufficient time for the transition to state 2. This velocity is 
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always larger than the particle velocity vq, as vq/vh = t^i/ii^i + ^,2) < 1- Thus the current, 
which is given by J = K2 ( 1 — p) for p ?a 1 , drops more steeply for p 1 than it increases for 
small p. 

While the current is the most important characteristics of a traffic system from a the- 
oretical point of view, it is not easily accessible in experiments. A quantity that can be 
measured directly is the motor velocity, which can be determined by tracking individual 
labeled motors in a background of unlabeled motors [2^. The velocity v can be expressed 
in terms of the current J using the relation v = Ji/ p. The density-dependence of the motor 
velocity is shown in Fig. [2](b). In contrast to the usual ASEP, which is obtained in the limit 
Ki ^ K,2 and for which the velocity exhibits a linear decrease as a function of the density 
p, the velocity is a convex function of the density in the two-step model. The reduction of 
the motor velocity is thus less pronounced at low motor densities, in particular if the actual 
step is not rate-limiting, i.e. if k,2 ^ ^^i- This means that a strong slowing down of the 
motors due to the traffic may only be observed for very dense motor traffic. At intermediate 
densities, the effect of exclusion is reduced for large k,2- This can be interpreted as a par- 
tially synchronized motion of motors: If one or several motors in state 2 are directly behind 
a motor in state 1, they follow almost immediately when the leading motor steps forward. 

The density- dependence of the motor current and the motor velocity can be described 
analytically using a mean-field approximation: all traffic effects can be subsumed into an 
effective rate for the actual stepping by replacing the rate k,2 by /t2(l — p). Then the 
probabilities Pi and P2 that the motor is in conformation 1 or 2 are calculated as for a 
single motor. This approximation leads to 

Pi= ^'^^-P^ , and P2 = ^ (1) 

Ki + K2[l-p) Ki + K2[l-p) 



and thus to 



and 



J p = PK2P2 = : 7z r (2) 

Kl + K2{1 - p) 



Vip) = n2iP2 = 

Ki + K2{1 - p) 

These expressions are plotted in Fig. Wis) ^^'^ i'^)-! again normalized by the single motor 
velocity vq. These plots show qualitatively the same behavior as the simulation data in 
Fig. [2](a) and (b). However, except for the limiting case ki ^ ^2-, which corresponds to the 
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usual ASEP, the mean field approximation overestimates the effects of the conformational 
transition. In particular, in the limit k.2 ^ /^i it predicts a linear increase of the current, 
J = {vq/£)p, and a constant velocity, f = t^o; up to almost the maximal density p = 1. The 
maximal current thus approaches Ji/vQ = 1 and occurs at a density value that approaches 
p = 1. In contrast, the density- dependence of the current obtained from simulations shows 
pronounced deviations from the linear density dependence of the current even for low ki/k,2 
and appears to converge to a function with a maximum of J/ (vo/i) ~ 0.37 at p ^ 0.75. The 
curves for ki/k,2 = 0.1 in Fig. [2] are close to this limit. 



III. MODEL WITH BINDING/UNBINDING KINETICS 

Next, we consider the effect of motor binding to and unbinding from the filament 
[Fig. [Hb)]. Even the simplest model with only two internal states of the motor allows 
many different implementations of binding and unbinding: unbinding may occur from state 
1 and/or state 2, described by the rates ei and €2, and unbound motors may also become 
bound in either state, described by the rates nic and 7i2C, where tti 2 are the second-order 
binding rates and c is the concentration of unbound motors. For simplicity, we restrict the 
discussion to the two cases, where unbinding occurs predominantly in one of the two states. 
Furthermore, our simulations indicate that whether newly bound motors are in state 1 or 
state 2 has only a small effect as l ong as the motors are processive, i.e. if they typically 



make many steps before unbinding 25 1. 



Results from simulations that include binding and unbinding of motors are shown in 
Fig. [31 The solid and dashed curves in Fig. [3] correspond to cases where unbinding occurs 



predominantly from state 1, with ei = Ki/100 and 62 = ei/100 26|], and only from state 2, 
with €2 = K2/IOO and ei = 0, respectively. In both cases, the unbinding rates are chosen to 
match a run length of 100 steps for the single motor. For processive motors, including motor 
binding and unbinding has only a small effect on the density-dependence of the motor current 
[Fig. [3] (a)] and the motor velocity compared to the case without binding/unbinding (dotted 
line). The density p of bound motors, however, is now determined by the binding/unbinding 
equilibrium. Fig. [3^b) shows the bound motor density as a function of the unbound motor 
concentration c, which is normalized by the desorption constant K = €0/(711 + 712). Here 
eo is the effective unbinding rate of a single motor, which is obtained from the general 
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effective unbinding rate ecs = ei-Pi + e2-P2 by taking the limit p = 0, which leads to eo = 
(eiK2+e2/«i)/(Ki + K2)- These curves are compared to simple Langmuir-type binding with the 
same effective unbinding rate eo for a single motor (dotted line), as used in models without 
internal states. While all three curves are the same for small binding rates (or unbound 
motor concentrations), there are marked differences for high binding rates. If unbinding 
occurs from state 1, binding is clearly stronger than in simple Langmuir-type binding 
while it is weaker if unbinding occurs from state 2. This behavior is not surprising, because, 
in denser traffic, motors spend a larger fraction of time in state 2, which increases the 
effective unbinding rate ees if unbinding occurs from state 2, but decreases it if unbinding 
occurs from state 1, as shown in Fig. [31(c). 

Closely connected with the unbinding rate is the run length or processivity, i.e. the num- 
ber of steps a motor takes before detaching from the filament, which is shown in Fig. ^d) 
as a function of the motor density on the filament. For simple Poissonian steppers, the 
unbinding rate is independent of the motor density, but the run length decreases linearly 
with increasing density due to the linear decrease of the velocity. For motors with a con- 
formational stepping cycle, the run length also decreases if unbinding occurs from state 2 
[although the decrease is not necessarily linear, in particular if ki ^ ^2, due to the con- 
vex density-dependence of the velocity as shown in Fig. EJ^b)]. If unbinding occurs however 
predominantly from state 1, the run length remains unaffected by the density of the motor 
traffic up to very high densities (where the small unbinding rate from state 2 becomes dom- 
inant). In this case, unbinding is characterized by a constant unbinding probability per step 
rather than a constant unbinding rate as in the case of Poissonian steppers. 



IV. EXPERIMENTAL RELEVANCE 



We have discussed the effects of the chemomechanical stepping cycle of a molecular motor 
on the traffic of many molecular motors using a simple two-state chemical cycle. For this 
model we have determined observable quantities, in particular the motor velocity and the 
average run length, which can be measured by tracking individual motors on a crowded 
filament |24|. 

Our analysis shows that the more complex stepping kinetics tend to diminish the slowing 
down of motors due to traffic congestion, but this effect is not very large. It is possible 
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that this effect contributes to the surprisingly small decrease of the velocity observed in the 



experiment of Seitz and Surrey 



2J]. In this experiment, the velocity of kinesin 1 motors 



remained almost constant up to unbound motor concentrations that resulted in nominal 
bound motor densities of ~0.3-0.5 and reduced the motor binding rate about two-fold. The 
interpretation of this experiment is however complicated by another observation showing 
that microtubules saturated with motors exhibit an approximately two-fold higher motor 
density in the presence of a non-hydrolysable ATP analogous than in the presence of ATP. It 
is thus rather likely that the specific geometry of binding of motors to the filament also plays 
an important role, in particular whether they bind with one or two heads and how many 
sites a motor occupies on the filament 2^. A larger size of the particles also diminishes the 



14, 



decrease of the velocity 



281 ]. so one may speculate that the relatively constant velocity 
observed by Seitz and Surrey is due to a combination of (at least) two weak effects. 

A more pronounced effect of the cycle is seen in the density-dependence of the run length 
and the effective unbinding rate. Depending on the state from which unbinding from the 
filament predominantly occurs, the effective unbinding rate may both increase or decrease 
with increasing traffic density, while the run length decreases or remains essentially constant. 



The latter behavior has been observed in the experiment of ref. [2J]. In this experiment, 
the time a motor spends bound to the filament before unbinding was found to increase 
about two-fold over the studied range of motor densities, indicating that the unbinding rate 
decreases in dense traffic, while the run length remained approximately constant. Another 
study [29| however reported the exact opposite, namely that unbinding is increased due to 
crowding. The reason for the discrepancy is not known, but there are several important 
differences in the experimental conditions. In the experiment of ref. 29| , microtubules were 
covered by very high densities of kinesin (up to one kinesin dimer per microtubule binding 
site) in the absence of nucleotides and unbinding was measured indirectly by light scattering 
during the transient dynamics after the addition of ATP. In contrast, the experiment of ref. 
2J] used lower motor densities, addressed the steady-state of the motor traffic, and motor 
properties were measured directly by tracking individual motors. The observation that 
different amounts of kinesin are needed to saturate microtubules in the presence of different 
nucleotides (ATP versus non-hydrolysable analogue) 2J] also suggests that kinesins may be 
bound in different ways in the two experiments. In general, even the equilibrium binding 
of dimeric motors to filaments (in the absence of ATP, i.e. without active movements) may 
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exhibit a range of stoichiometries and rather complex dynamics 30 1. 

Finally, the main effect of crowded molecular motor traffic is the emergence of a traffic- 



jam-like domain of 
several experiments 



ligh motor density at the end of a filament 



13 



21 



231 ] ■ as observed in 



22| . If the filament is in contact with a large reservoir of unbound 



motors, the length of the high density domain or traffic jam is of the order of the run length 
in ASEP-type models 2|. If the unbinding rate is reduced in dense traffic, as in our model 



with unbinding from state 1 and as suggested by the experiment of ref. 



24l |. one may expect 



to observe much longer traffic jams. Extending the arguments given in ref. Q|, one obtains 
the estimate that the jam length is larger than the run length by a factor eo/e(p). 



Appendix 

The model we have discussed in this article represents the simplest possible chemome- 
chanical cycle of a molecular motor with two subsequent transitions that are treated as 
irreversible. This minimal model of a stepping cycle provides a good description of the dy- 
namics of kinesin motors under typical experimental conditions (low concentrations of the 



hydrolysis products ADP and Pi, no opposing force), see, e.g., Ref. [2^, where backward 



steps are rare 



18 



3l| , but should not be expected to be valid under all possible experimental 



conditions or for other types of motors such as dyneins which exhibit backward steps more 



frequently than kinesins 32l, l33|] . On the other hand, it is straightforward to include back- 
ward steps in the theoretical framework used here as already described in Ref. jsl. To test 
to what extent our results remain valid if the dynamics of motor stepping are more complex, 
we simulated several variants of our model that include reverse transitions and backward 
steps. All these variants exhibit qualitatively the same behavior as the simple model. 

Surprisingly, the main difference to the simple model is already obtained when the reverse 
transition for the transition from state 1 to state 2 is included: In the original model with 
irreversible transitions, a motor in dense traffic often waits in state 2 until the site in front 
becomes available. If the reverse transition is included such a motor will switch between 
states 1 and 2 while waiting, so the probability that the motor is in state 2 is reduced in dense 
traffic compared to the irreversible model. As a consequence the motor velocity is slightly 
reduced and the distinction between the cases with unbinding from state 1 or 2 becomes less 
pronounced if the reverse rate is sufficiently large. In particular, the run length exhibits a 
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weak decrease with increasing motor density (rather than being constant) in the case where 
unbinding occurs mainly from state 1. Essentially the same behavior was observed for a 
cycle that included reverse transition for both the internal transition and the actual step. 

Finally we also simulated the case where a motor in state 2 may make either a forward 
step or (with a smaller rate) a backward step, which provides the simplest model where 
forward and backward steps occur along different reaction pathways, a situation suggested 
for kinesin by both experiments and modeling [sil . Our simulation results for this case 
strongly resemble those for the simplest cycle without backsteps or reverse transitions. 
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FIG. 1: Model for traffic of molecular motors with a two-state chemomechanical cycle: (a) Motors 
can be in two conformational states, state 1 (white) and 2 (grey). The two states may be considered 
as a weakly bound state and a strongly bound state. A motor in state 1 can make a transition to 
state 2 with rate ki, while remaining at the same site on the filament. A motor in state 2 makes 
a transition to state 1 with rate K2. The latter transition is accompanied by a step of size ^ to the 
next site on the filament and occurs only if the forward neighbor site is not occupied by another 
motor in either state 1 or 2. (b) Motors in state 1 and 2 unbind from the filament with rates ei 
and £2, respectively. Likewise, a free binding site on the filament becomes occupied by a motor in 
state 1 or 2 with rates vric and 7r2C, respectively, which depend on the concentration c of unbound 
motors. Unbound motors may not bind to sites already occupied by a motor in either state. 
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FIG. 2: Effect of the chemomechanical stepping cycle: (a,c) Motor current J and (b,d) velocity 
V as functions of the motor density p on the filament for different choices of the rates ki and K2- 
(a) and (b) show results from simulations, while (c) and (d) show the corresponding mean field 
results. All plots shown here have been obtained for the simplified models in which the motors do 
not unbind from the filament. Both current and velocity have been normalized by the single motor 
velocity vq, which is taken to be an experimentally determined quantity. With this normalization 
these quantities are functions of the ratio of 1^1/1^2 only. 
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FIG. 3: Effect of motor binding and unbinding: (a) Motor current J as a function of the motor 
density /), (b) motor density /? as a function of the unbound motor concentration c, normahzed 
by the dissociation constant K = eo/^ni + tt2), (c) effective unbinding rate egfj and (d) average 
run length Axb as functions of the bound motor density p. All plots show simulation results for 
models with unbinding predominantly from state 1 (solid lines, ei = ki/100, e2 = ei/100, ki = K2) 
and with unbinding from state 2 (dashed lines, ei = 0, £2 = K2/IOO, ki = K2). The dotted 
line in (a) shows the corresponding curve without binding/unbinding, the dotted line in (b) the 
Langmuir-type binding curve for the same single motor unbinding rate cq. 
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